Temporal decorrelations in compressible isotropic turbulence 
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Temporal decorrelations in compressible isotropic turbulence are studied using the space-time 
correlation theory and direct numerical simulation. A swept-wave model is developed for dilatational 
components while the classic random sweeping model is proposed for solenoidal components. The 
swept-wave model shows that the temporal decorrelations in dilatational fluctuations are dominated 
by two physical processes: random sweeping and wave propagation. These models are supported 
by the direct numerical simulation of compressible isotropic turbulence, in the sense of that all 
curves of normalized time correlations for different wavenumbers collapse into a single one using the 
normalized time separations. The swept-wave model is further extended to account for a constant 
mean velocity. 

PACS numbers: 47.27.eb, 47.27.Gs, 47.40.-x 
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A milestone in isotropic and homogeneous turbulence 
is the random sweeping hypothesis P, 0- The random 
sweeping hypothesis proposes that there is a temporal 
decorrelation process in incompressible isotropic turbu- 
lence and develops a simple model for space-time cor- 
relations of velocity fluctuations These results are 
examined theoretically 3-5] and verified experimentally 
Q and numerically [3, Hj]. The space-time correlation 
models are used to predict the scalings of wavenumber 
or frequency energy spectra in turbulent flows [j| [lfjj |. 
The decorrelation processes are also relevant to the non- 
Gaussian statistics [llj and intermittency Their 
further applications can be found in turbulence gener- 
ated noise 13J. The recently increasing studies on com- 
pressible isotropic turbulence raise such a question on 
the effects of compressibility on decorrelation processes 
1J-|16|. In this letter, we will study the decorrelation 
processes in compressible isotropic turbulence and pro- 
pose a model for space-time correlations of dilatational 
components. 

A compressible turbulence is associated with two char- 
acteristic velocities: fluid velocity and sound speed, 
whereas an incompressible one is only associated with 
fluid velocity. Therefore, the decorrelation processes in 
compressible turbulence are very different from incom- 
pressible one. A space-time correlation is the essential 
quantity to measure the decorrelation processes in turbu- 
lent flows. Three typical models exist for space-time cor- 
relations in turbulence theory. The first one is, as stated 
above, the random sweeping model for incompressible 
turbulence We will show that it cannot character- 
ize the acoustic components in compressible turbulence. 
The second one is the Taylor frozen flow model 17|. It 



has been shown that this model is not a good approxi- 
mation for dilatational components [18fl. The third one 
is the linear wave propagation model [18|. This model 
can be used for dilatational components if compressible 
turbulence has a dominating mean velocity. However, we 
will find that it does not decrease with increasing tem- 
poral separation, which violates the nature of correlation 



functions. 

In the present letter, we will develop a space-time 
correlation model for compressible isotropic turbulence. 
This is achieved by the Helmholtz decomposition: a ve- 
locity field can be split into the solenoidal and dilata- 
tional components. A swept-wave model will be devel- 
oped for the dilatational components while the solenoidal 
components are expected to follow the random sweeping 
model. The swept-wave model will be numerically val- 
idated and further used to elucidate the decorrelation 
process in compressible turbulence. 

We consider compressible and isotropic turbulence 
with periodic boundary conditions. In this case, the 
Helmholtz decomposition for velocity fields can be made 
as follows 10] 
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where u s and u d are the solenoidal (i.e. incompress- 
ible) and dilatational components, respectively. The har- 
monic component is taken to be zero. We will investigate 
the temporal decorrelations of solenoidal and dilatational 
components. The temporal decorrelations can be mea- 
sured by the space-time correlations of velocity fluctua- 
tions 



R(r,r) = (u e (x,t)ui(x + r.t + T)), 
or its equivalent forms in Fourier space 

R(k,r) = (u e (\s.,t)u e (-k,t + T)). 



(2) 



(3) 



Here, r and k are the magnitudes of separation vector 
r and wavenumber vector k. The similar quantities R s 
and R d can be defined for the solenoidal and dilatational 
components u s and u d , respectively. 

We propose that a solenoidal component follows the 
same decorrelation process as the random sweeping pro- 
cess for incompressible isotropic turbulence: small eddies 
are randomly convected or swept by energy-contained ed- 
dies, where the contribution of dilatational components 
to the energy-contained eddies is comparably small [3| . 
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FIG. 1. Normalized time correlations of solenoidal compo- 
nents vs time separations for k = 10, 20, 30, 40, 50, 60. 



The random sweeping process can be described by a sim- 
ple idealized convection equation as follows [l[ 



(4) 



where v = (vi, V2, V3) is a spatially uniform and station- 
ary Gaussian random field. V — |v|/\/3 is the r.m.s of 
solenoidal components along one axis. Its solution in the 
Fourier space is given by 



fif (k, t) = uj(k, 0) exp [-i(k • v)i] . (5) 
The time correlation of Fourier mode is formulated as 



R s (k, t) = R s (k, 0) exp — V 2 k 2 r 2 



1. 



(6) 



A dilatational component in compressible isotropic 
turbulence propagates at the speed of sound relative to 
moving fluids. This implies that the dilatational fluctu- 
ations are swept by the energy-contained eddies. There- 
fore, the temporal decorrelations in dilatational compo- 
nents are governed by two dynamic processes: random 
sweeping and wave propagation. The well-known lin- 
ear wave propagation model 18] only includes the wave 
propagation process. In order to account for the random 
sweeping effect, we introduce a new term vV into the 
linear wave propagation equation and propose the gov- 
erning equation for dilatational fluctuations as follows 



D 2 



- a 2 V 2 )uf = 0, 



(7) 



FIG. 2. Normalized time correlations of solenoidal com- 
ponents vs normalized time separations Vkr for k — 
10, 20, 30, 40, 50, 60. 



where a is the mean speed of sound and 

d fa „\ fa a . 
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Here, v is the same Gaussian random field as one in 
Eq. (jlj). The new term vV in Eq. (JJJ represents the ran- 
dom sweeping effect, which is absent in the linear wave 
propagation model [18| ■ The solution of Eq. in Fourier 
space is given by 

u|(k, t) = u^ + (k, 0) exp H(k • v)t - ikdt] + 

uf~(k,0)exp[-i(k.v)t + ifcat], (9) 

where and ii d t ~ are the Fourier coefficients. 

The time correlation of Fourier mode is calculated as 
follows 



R d (k,r) = (uf(k,t)uf(-k,t + T)) 



(10) 



(k, 0)u^ + (-k, 0))(exp[z(k • v)r + ikar}} 
+ (u d ~ (k, 0)uf- (-k, Q))(exp[« (k • v)r - ikar}). 

Here, the mode correlation R d (k, 0) is given by 

l -R d (k,Q) = (u d+ (-k,0)u d+ (k,0)) = (fi|-(-k,0)^-(k,0)). 

Therefore, the correlation functions can be expressed as 



R d {k, T ) 
R d (k,0) 



cos (fear) exp 



k 2 V 2 r 2 



(11) 



The swept-wave model (fTTj) contains two factors: a lin- 
ear wave function and an exponential function. The 
first factor represents the random sweeping effect and the 
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second one represents the wave propagation process. If 
the sweeping velocity is zero, it becomes the linear wave 
propagation model. In fact, the linear wave propagation 
model in compressible isotropic turbulence is simplified 
as cos(fcar), which is the inverse Fourier transformation 
of Equation (19) in [lij . This cosine function does not 
decay to zero as time separation increases. 

To validate the swept-wave model, we solve the three- 
dimensional, compressible Navier-Stokes equations in a 
cubic box of side 2n with specific heat ratio 7 = 1.4 and 
Prandtl number Pr — 0.7 using an optimized sixth order 
compact, finite difference scheme. Statistically station- 
ary flow fields are achieved by including a forcing term 
f(x, t) with solenoidal modes only, V • f = 0. The forc- 
ing term is nonzero only in the range 1 < k < 3 and 
obeys a Gaussian random distribution with an exponen- 
tial temporal correlation. Each component of f(k,i) is 
defined by fi = (Sij — kikj/k 2 )cjj, where jjj is generated 
by an independent Uhlenbeck-Ornstien process. After 
the flow fields become statistically stationary, a total of 
400 flow fields with computational time increment 0.02 
are chosen to calculate time correlations. The Taylor's 
micro-scale based Reynolds number is about 80 and the 
turbulent Mach number is about 0.63. We also compare 
the space-time correlations from the present case with 
ones from decaying turbulence. The results obtained are 
consistent. 

Fig. [1] shows the correlation coefficients of solenoidal 
modes for wavenumbers k = 10, 20, • • • 50 and 60, span- 
ning a range of scales from the integral scale to the 
dissipation scales. The correlation coefficients are the 
normalized correlation functions by the mode correla- 
tion R s (k,0). Obviously, the solenoidal modes decorre- 
late more quickly at larger wavenumbers than at small 
wavenumbers. These results in Fig. [1] are all plotted to- 
gether in Fig. [2l with the horizontal axis defined by the 
normalized time scale Vkr. This normalization causes 
excellent collapse of the correlation coefficients. The col- 
lapse on the normalized time scale Vkr supports the ran- 
dom sweeping hypothesis for the solenoidal components. 

Fig. [3] plots the normalized time correlations 
R (k, t)/R (k, 0) for dilatational components from DNS 
data for wavenumber k — 10, 20, • • • and 60, where the 
correlations are normalized by the correlation R d (k,0). 
It is observed from Fig. [3] that the time correlations of 
dilatational components decay with oscillations. This is 
very different from solenoidal components where the time 
correlations decay without any oscillation. These oscil- 
latory decays confirm that temporal decorrelations in di- 
latational components are mainly determined by both 
random sweeping and wave propagation. 

Fig. [4] presents the normalized time correlations 
R (k, t)/R (k, 0) versus the normalized time separation 
Vkr for wavenumber k = 10, 20, • • • , 50 and 60. The time 
normalization leads to the virtual collapse of all curves. 
The collapsed curves decay with oscillations. This re- 
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FIG. 3. Normalized time correlations of dilatational compo- 
nents vs time separations for k — 10, 20, 30, 40, 50, 60. 
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FIG. 4. Normalized time correlations of dilatational com- 
ponents vs normalized time separations Vkr for k = 
10, 20, 30, 40, 50, 60. 



suit verifies the proposed swept-wave model for dilata- 
tional components. We note that the ratio of two scale- 
similarity variables Vkr and akr is constant. 

Fig. [5] compares the collapsed curves for solenoidal 
components with the collapsed ones for dilatational com- 
ponents, where the horizontal axis is normalized as the 
scale-dependent similarity variable Vkr. It is observed 
that the collapsed curves for solenoidal components act 
as an envelop of the collapsed one for dilatational compo- 
nents. This confirms the swept-wave model where the ex- 
ponential function exp(— 0.5V r2 fc 2 r 2 ) acts as an envelop. 
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FIG. 5. Normalized time correlations vs normalized 
time separations Vkr. Dilatational components for k — 
10, 20, 30, 40, 50, 60; solenoidal component and its mirror im- 
age for k = 30. 



propagation while a solenoidal component dominated by 
the random sweeping effect. We further develop a swept- 
wave model for dilatational fluctuations. This model is 
distinct from the linear wave propagation model since 
it includes the random sweeping process. The DNS data 
validates the swept- wave model for compressible isotropic 
turbulence. The further extension of the swept- wave 
model to turbulent shear flows is referred to the elliptic 
model 2l| and the present model can be used to study 
the propagation velocity of coherent structures in com- 
pressible turbulence. 
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We can further calculate the space-time correlations of 
dilatational components using the Fourier transformation 
of R d {k, t) from wavenumber space to spatial one 



V 2 k 2 r 2 



_ sin(fcr) 

cos(fear) — dk. 

kr 



(12) 



Eq. (1121) can be extended to account for a constant mean 
velocity. Without loss of generality, we choose the con- 
stant mean velocity U\ in the direction of the x\ axis. Ap- 
plying the coordinate transformation (yi = xi~U\t,y2 = 
£2,2/3 = Z3) to Eq. (JTJ), we obtain 



R d (r, t) = E(k)exp (-\v 2 k 



2 ;„2 T 2 



cos(fcar) S ^ n , ^ r T~~ dk. 
k{r — U\t) 



(13) 



In comparison with the linear wave propagation model, 
this model (|13[) contains an additional exponential func- 
tion that is responsible for the random sweeping effect. 
It also confirms that Taylor's frozen flow model is not a 
good approximation to the space-time correlation of di- 
latational components. Wilczek and Narita [2(| consider 
the random sweeping model with constant mean velocity. 
The present model is consistent with their results. 

In summary, we find that solenoidal and dilatational 
components in compressible isotropic turbulence display 
different decorrelation processes: a dilatational compo- 
nent is dominated by both random sweeping and wave 
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